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ABSTRACT: In this paper we have studied die motion dad stability of the centre of mass of a system of two satellites 
connected by inextensible calve under the influence of air resistance and magnetic force in the central gravitational field of 
oblate earth in circular orbit. W'c have obtained an equilibrium point which has been .shown in he stable in the sense of 
Liapunov. 
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I. INTRODUCTION 

This paper is devoted to study the equilibrium position under the influence of air resistance and magnetic force of 
oblate earth in case of circular orbit of the centre of mass of the system. 

For this, firstly we have derived equations of motion in case of circular orbit of the centre of mass of the system 
under perturbative forces mentioned above and then Jacobian integral for the problem is obtained. Equilibrium Position has 
been obtained shown to be stable in the sense of Liapunov. This work is direct generalization of works done by V.V. 
Beletsky; R. B. Singh; B. Sharma; S. K. Das; P. K. Bhattacharyya and C.P.Singh. 

II. EQUATIONS OF MOTION OF ONE OF THE TWO SATELLITES IN ELLIPTIC ORBIT 

The equations of motion of one of the two satellites when the centre of mass moves along a keplerian elliptical orbit 
in Nechvill's co-ordinate system have been derived in the form. 
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Where, 
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/toe = A x = : A being Lagrange's multiplier's and mi, m 2 being masses of two satellites. 

/u jU l Q niyn 2 

l Q being the length of cable connected by two satellites 
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p = — = ; p being focal parameter and e eccentricity of the orbit of centre of mass 

p 1 + ecosv 

R = Radius vector of the centre of mass from the attracting centre 
v = True anomaly of the centre of mass 

i = Inclination of the orbit of centre of mass with the equatorial plane of the earth 

A = — — = oblateness force parameter 
P 
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B = -. = magnetic force parameter 
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f = . = Air resistance force parameter 
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Here, dashes denote differentiations w.r. to true anomaly v. 
The condition of constraint is given by 

x 2 + y 2 + z z < (2.2) 

p[ 
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Since the general solution of the system of differential equations given by (2.1) is beyond our reach so we restrict ourselves 
to the case of circular orbit of the centre of mass of the system in equatorial plane (i = 0) and hence we get from (2.1) on 



putting. 

p = = 1 , p ' = 0 and i = 0 (for equatorial plane) 

1 + ecosv 
x"~ 2y '- 3x = I a x - 4Ax - B 
y "+ 2x ' = X a y + Ay - f 

z"+Z = AaZ + Az (2.3) 

The condition of constraint given by (2.2) takes the form 

x 2 + y 2 + z 2 < 1 (2.4) 

Thus if inequality sign holds in (2.4) then the motion takes place with loose string and the motion is called free motion. If the 
equality sign holds in (2.4), then the motion takes place with tight string and the motion is called constrained motion. 
We are actually interested in stability due to constrained motion. Thus, motion takes place on unit sphere given by - 

x 2 + y 2 + z 2 =l (2.5) 

Differentiating (2.5), we get 

xx l + yy l + zz = 0 (2.6) 

To obtain Jacobian integral of the problem, we multiply the first, second and third equation of (2.3) by x', y' and z' 
respectively and add them together, we get after integrating on using (2.5) and (2.6) 

x ' 2 + y a + z a - (3x 2 - z 2 ) = 5 Ax 2 -2Bx-2Jy + h (2.7) 

Which is known as Jacobian integral and can be interpreted as energy equation with modified potential given by 

V J(3r'-- , )-^.v-/>'.v-/v (2.8) 

Differentiating (6) with respect to v, we get 

x a +y' 2 +z a =-(xx"+yy"+zz") (2.9) 

Multiplying the first, second and the third equations of (2.3) by x, y and z respectively and adding we get on using (2.5) 

xx"+ yy"+ zz" = 2(xy'- x' y) + (3x 2 - z 2 )-(5x 2 -1)a+J x - Bx- fy (2.10) 

Using (2.9) in (2.10), we get 

-Aoc = (x' 2 +y' 2 +z a ) + 2(xy'-x'y) + (5x 2 -l) + (3x 2 -z 2 )-Bx-jy. (2.11) 

To simplify (2.7) and (2.8), we use spherical polar coordinate on unit sphere: 

x = cos<f>cosy/:y = cos^sin^- andz = sin^ (2.12) 



Using (2.12), (2.7) and (2.8) become respectively 

(f) a +y/ a cos 2 (j) = ^3 cos 2 y/ + l^cos 2 ^ — 5Acos 2 ^cos 2 y/ -2Bcos<f>cosy/ -2f cos (f) sin yz + J^ 

(2.13) 

Where, 

\ = h - 1 and 

V(y,^) = -(3cos 2 ^/- + l)cos 2 ^ + 5Acos 2 ^cos 2 ^ + 25 cos ^ cos ^ + 2/ cos ^ sin (2.14) 



III. PARTICULAR SOLUTION AND STABILITY 

For equilibrium positions, (2.14) can be taken as modified potential energy. 

The equilibrium positions are given by the stationary values of v(j), y/^ and hence, we have 

dy/ 

Differentiating (2.14) partially with respect to (j) and using (3.1), we get 
Sin^ =0i.e. (j) =0 
Also, we have 
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= 3cos^ 0 sin^ () -5Acos^ 0 sin^ 0 -2?sin^ + /cos^ =0 

For smallest value of lf/ 0 , we have from (3.4) on putting COS^ 0 =1 and sin y/ Q = y/ Q 

3-5A-B 
Thus, the equilibrium point is given by 



$ = = 0 



-f 



3-5A-B 

To test the stability to the equilibrium position given by (2.6), we have 

\ d 2 V 

w 



= 2^3cos y/ Q + lJ-10Acos y/ - 2B cos y/ 0 - 2f sin y 
= 6 cos 2y/ Q -10 A cos 2y 0 - 2B cos y/ 0 -2f sin y 0 . 



e 2 v 



= 0 = 



5 2 V 



For equilibrium point given by (3.6) to be stable in the sense of Liapunov, we have to show that 

Using (3.7), (3.8) and (3.9) in (3.10) it can be easily seen that (3.10) is positive if 5A + B < 3. . 

Conclusion: we conclude that the equilibrium point 

<A = d () = Q; w = w = 1- 

ro y ro 3 _ 5A _ B 

is stable in the sense of Liapunov if 5 A + B < 3. 
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